Electronic Correlation effects in superconducting picene from ab-initio calculations 
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We show, by means of ab-initio calculations, that electron-electron correlations play an impor- 
tant role in potassium-doped picene (Jf^-picene), recently characterized as a superconductor with 
T c = UK. The inclusion of exchange interactions by means of hybrid functionals reproduces the 
correct gap for the undoped compound and predicts an antiferromagnetic state for x — 3, where 
superconductivity has been observed. These calculations, which do not require to assume a value 
for the interaction strength, indirectly suggest that these materials should have a sizeable ratio be- 
tween the effective Coulomb repulsion U and the bandwidth. This is fully compatible with simple 
estimates of this ratio. Using these values of U in a simple effective Hubbard model, an antiferro- 
magnetic state is indeed stabilized. Our results highlight the similarity between potassium-doped 
picene and alkali-doped fulleride superconductors. 

PACS numbers: 71.20.Tx,74.20.Pq, 74.70. Wz,74.70.Kn 
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I. INTRODUCTION 

The unanticipated discovery of superconductivity with 
critical temperature T c — 18 A" in an aromatic compound 
like potassium-doped picene^ opens a new path to super- 
conductivity in organic materials. Looking for an answer 
to the most natural question, namely "What is the origin 
of the electron pairing?" , we are tempted to consider only 
two mutually exclusive options: On one side standard 
electron-phonon superconductors which are expected to 
have low critical temperatures, on the other "anomalous" 
superconductors, dominated by electron-electron correla- 
tion effects. In this latter case pairing is assumed to have 
electronic origin while the critical temperature may even 
exceed 100K. 

A challenge to the above distinction comes from a fam- 
ily of organic superconductors, the alkali-metal-doped 
fullerides. There is wide consensus that pairing is me- 
diated by intramolecular phonons in these materials 2 . 
Nonetheless CS3C60, which has one of the largest in- 
tramolecular distances in the family, is an antiferro- 
magnetic insulator which turns into a 38K supercon- 
ductor under pressure 3 . This evidence confirms a pre- 
vious proposal that correlations are indeed important 
and helpful for superconductivity in fullerides^— despite 
the phononic nature of the pairing " glue" . In this con- 
text the discovery of relatively-high T c in a different or- 
ganic compound doped with potassium is exciting and 
brings, among many others, the question about the role 
of electron-electron correlations. This paper is devoted 
to a first investigation of the fingerprints of electronic 
correlations in i^ x -picene using a wide range of density- 
functional theory (DFT) methods, including in particu- 
lar hybrid functionals designed to include some degree of 
exchange correlations. 

Picene is an aromatic molecule (C22H14) composed by 
five benzene rings arranged in an "armchair" structure 



which determines a high chemical stability reflected in 
the band gap of 3.3 eV. Superconductivity has been dis- 
covered in solid picene with a monoclinic structure when 
the band insulating stoichiometric compound is doped 
by alkali atoms. Indicating with K K -picene a doped solid 
with x potassium atoms per picene molecule, supercon- 
ductivity has been observed in the region 2.6 < x < 3.3^. 
The first ab-initio study of K x -picene has been presented 
mi, where the electronic structure and Fermi surface 
have been obtained within the Local Density Approxi- 
mation (LDA) . In this work we extend this analysis using 
DFT-based methods aiming to include correlation effects. 
Namely we compare Perdew-Burke-Ernzcrhof (PBE) 8 
implementation of generalized gradient corrections with 
LDA+U 9 calculations for a tight-binding model derived 
from PBE, and we use Heyd-Scuseria-Ernzerhof (HSE) 10 
and B3LYP 11 hybrid functionals including different de- 
grees of many-body exchange. 

The paper is organized as follows: Section II presents 
results for undoped solid picene, while in Section III we 
turn to describe the doped solid picene estimating a siz- 
able on-site electron-electron interaction, mapping the 
doped solid picene into an Hubbard model, and then dis- 
cussing the emergence of an AFM solution. 



II. BANDSTRUCTURE OF UNDOPED PICENE 
A. Single Picene Molecule 

Given the molecular nature of solid picene, it is use- 
ful to recall the basis aspects of the molecular electronic 
structure, which is expected is the main building block 
of the bandstructure of the solid. Calculations for an iso- 
lated molecule are based on the GAMESS package^ and 
use PBE, PBEO and B3LYP functionals choosing TZV 
basis set. A PBE calculation fixing the molecule in its 
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equilibrium geometry in the solid gives a HOMO-LUMO 
gap of 2.9 eV as reported ini. The gap increases to 4.2 
(4.4) eV within B3LYP (PBEO), suggesting an impor- 
tant role of intramolecular exchange interactions. Since 
our focus is on the electron-doped system, the structure 
of the unoccupied states is of particular importance, and 
it is characterized by two closely lying states (LUMO 
and LUMO+1). The difference between these levels is 
0.19 (0.20) eV using PBE (B3LYP) and it can shrink 
down to 75 meV if the geometry of the molecule is op- 
timized, leading to a flatter structure which in turn re- 
duces the molecular overcrowding. On the other hand, 
using the geometry of the doped crystals, the LUMO and 
LUMO+1 gap increases to 0.36 eV. 

The potassium electrons will therefore populate the 
bands arising from the LUMO and LUMO+1, which can 
host up to four electrons. We briefly anticipate two possi- 
ble configurations for x — 3, the potassium concentration 
around which superconductivity has been observed: 

• (i) the bands originating from LUMO and 
LUMO+1 may remain distinct with two electrons 
per molecule filling the LUMO band and the re- 
maining electron going in the LUMO+1 band, 
which becomes half-filled. 

• (ii) The two bands are not well separated, and the 
three electrons partially populate both bands so 
that the system behaves more like a 3/4- filled two- 
band system. 

While the former scenario would lead essentially to a 
single band half-filled Hubbard model which is likely 
to order antiferrogmagnetically, the latter situation is 
less favorable for antiferromagnetic ordering and it can 
give rise to a competition between different phases in- 
volving magnetic and charge ordering, as shown by the 
phase diagrams of 3/4-filled TMTTF and TMTSF or- 
ganic compounds 13 . Given the small molecular splitting 
both situations are in principle possible. We will show 
that, at least in the present DFT calculations, the first 
scenario seems to be preferred, even if a more accurate 
treatment of electron-electron correlations may revive the 
second possibility. 



B. Undoped Solid Picene: Generalized Gradient 
Approximation and Hybrid Functionals 

We now move to pristine solid picene, that we first 
study via PBE and the projector augmented wave (PAW) 
method 1 using the Vienna ab-initio simulation package 
(VASP)i£. The valence pseudo-wave-functions were ex- 
panded in a plane-wave basis set with a cutoff energy of 
400 eV and all the integrations in the Brillouin zone are 
performed with a Gaussian broadening a = 0.02 using 
a sampling grid of 5x7x3 k-points along lattice vectors 
a,b,c of the unit cell. The lattice parameters are fixed at 
the experimental values^. Within the unit cell there are 
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FIG. 1. (Color online) Schematic arrangement of two inequiv- 
alent picene molecules (A and B) along the a,b planar molec- 
ular axes in a) pristine and b) doped molecular crystals. 



two inequivalent picene molecules (A and B) arranged in 
herringbone structure in the ab plane with stacking along 
the c direction (see Fig. [T|). 

The bandstructure obtained with PBE confirms the 
molecular nature of the solid because the relevant bands 
clearly correspond to the molecular orbitals described 
above. Our results are in good agreement with Rcf. 7 , 
where the relation between the molecular orbitals and 
the maximally localized Wannier orbitals of the solid 
has been discussed in details. Since there are two in- 
equivalent picene molecules, the HOMOs of the picene 
molecules hybridize giving rise to two bands just below 
the Fermi level, which are separated by a band gap of 
2.2eV from a manifold of four bands originating from 
the two LUMOs and two LUMO+ls. As customary, 
PBE underestimates the experimental value of the gapii. 
Both the valence and the conduction bands have a band- 
width of around 0.5 eV. The dispersion is enhanced along 
the a*,b* directions as expected from the layered crystal 
structure. For a similar molecular crystal, Van der Waals 
interactions may in principle play a role in determining 
the structure. The reasonable agreement between the 
experimental structure and the results of our theoretical 
relaxation suggests that this role is indeed not crucial. 

To overcome the limitations of PBE (and LDA) and its 
shortcoming in the determination of the gap, we repeated 
the same kind of calculations using HSE and B3LYP hy- 
brid functionals, with the same basis set and computa- 
tional details mentioned above. The two approaches lead 
to an enhanced gap of 3.0 and 3.25 eV respectively. The 
inclusion of some degree of short-range exchange interac- 
tions (within a Hartree-Fock treatment) is therefore es- 
sential to reproduce the experimental gap of 3.3 eV, sig- 
naling the important role of intramolecular interactions. 
An inspection to the density of states obtained via HSE 
and B3LYP clarifies the role of of the exchange interac- 
tions in the undoped system. The bands become signif- 
icantly narrower with respect to PBE as a consequence 
of a more localized nature of the carriers, as expected in 
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FIG. 2. (Color online) Results for pristine Picene: Panel (a) 
presents the PBE bandstructure; panel (b) shows the density 
of states calculated with PBE, HSE, and B3LYP functionals. 
The energies are plotted along lines in the Brillouin zone con- 
necting the points T = (0,0,0), X = (±,0,0), M = (±, §,0), 
L = (±, 0, ^). The zero of the energy is at the Fermi-level. 



a correlated system. Correspondingly, the wavefunctions 
at each site become even closer to the molecular orbitals, 
another signature of strong correlations. This effect is 
particularly strong in the unoccupied bands coming from 
LUMO and LUMO+1 that the doped electrons are ex- 
pected to populate. 



III. BANDSTRUCTURE OF ^-PICENE 
A. PBE bandstructure 

We finally turn to K x -picene. The experimental results 
already show that a rigid band picture does not hold, as 
shown by the difference in the measured lattice param- 
eters of undoped and K^.g-picene, which suggests that 
the dopants are not intercalated in the interlayer region 
and that they affect more deeply the bandstructure, as 
already discussed in ReL 1,7 . We built a unit cell with 
the lattice parameters measured for K^.g-picene and we 
inserted 3 potassium dopants per picene molecule in the 
ab plane as suggested in Ref^i (see Fig. Q]). We also 
repeated within PBE an optimization of the potassium 
positions starting from different configuration obtaining 
results essentially equivalent to LDA optimization^. The 
dopant orbitals are higher in energy and very weakly 
hybridized with the molecular states close to the Fermi 
level. Therefore their outer s electrons are donated to the 
picene bands. Indeed the four bands originating from the 
LUMO and LUMO+I of the two inequivalent molecules 
arc filled by the six electrons donated by potassium atoms 
(three K per picene molecule). Within PBE (see Fig. [3^,) 
the LUMO and LUMO+I bands are slightly separated 
so that the two low-lying bands are completely filled, 
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FIG. 3. (Color online) Results for I<3-Picene restricted to the 
energy region close to the Fermi level: (a) shows the PBE 
band-structure, while (b) presents the density of states cal- 
culated with PBE (non magnetic), compared with HSE both 
in nonmagnetic and AFM state (b). The energies are plot- 
ted along lines in the Brillouin zone connecting the points 
r = (0,0,0), X = (±, 0,0), A/ = (i,±,0), L = (i,0,i). The 
zero of energy is at the Fermi-level. 

while the two LUMO+1 bands host an average of one 
electron per site per band, or, in other words, they are 
half-filled. The bandwidth W is of around 0.8 eV for 
all the relevant bands. The above described electronic 
structure, with one electron per picene molecule in the 
highest partially occupied band, can give rise to magnetic 
ordering of the spins of the electrons in the LUMO+1 
bands. Within PBE we have stabilized both a ferromag- 
netic solution and an antiferromagnetic (AFM) solution 
in which A and B molecules have respectively parallel 
and opposite spins. Both solutions retain metallic char- 
acter even if they present a momentum of around 0.25 
fiB- Their energy is comparable with the non magnetic 
state within the accuracy of the calculation. 



B. Antiferromagnetic solution within mean- field 

In order to obtain a more reliable picture of the mag- 
netic behavior of the system we need to better account 
for electron-electron interactions, which are in princi- 
ple able to give rise to a Mott insulating state with or 
without magnetic ordering. A simple way is to include 
the molecular Coulomb repulsion measured by the Hub- 
bard parameter U. It is therefore important to esti- 
mate the value of U, which is essentially the energetic 
cost to doubly occupy a molecule including the screen- 
ing effects of the other bandsi^. We can compute U 
for a single molecule with three electrons from the en- 
ergies of the molecule charged with 2,3 and 4 electrons 
as U = E{A) - 2E(3) + E(2), where E(M) is the to- 
tal energy of a molecule charged with M extra electrons, 
and we will consider N — 3 in the following. Such an 
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estimate for an isolated molecule needs to be corrected 
in order to include the screening effects in the solid. A 
first estimate can be obtained by considering the effect 
of the polarization of a charged molecule^ placed inside 
a cavity of an homogeneous dielectric medium charac- 
terized by a dielectric constant e. Using typical values 
for organic molecular crystals e = 3,4, 5 and 6, we ob- 
tain respectively U = 1.67, 1.43, 1.29, 1.20 eVSS. Even if 
this estimate can not be regarded as quantitative, given 
the bandwidth of W ~ 0.8eV we obtain a large value of 
the ratio U/W, very close to similar estimates for doped 
Ceo, where the signatures of the strong repulsion have 
been experimentally observed 3 . It is therefore impor- 
tant to explore the consequences of electronic correlations 
also in picene-based solids. A first approach is to intro- 
duce the Hubbard repulsion at a mean-field level using 
the DFT bands as the bare bandstructure. Therefore we 
built a tight-binding representation of the bandstructure 
using Wannier90 21 to compute the maximally localized 
Wannier orbitals starting from the PBE bandstructure 
obtained using Quantum Espressc 22 . Within the tight- 
binding representation the system is represented as a lat- 
tice model, in which the lattice sites coincide with the 
picene molecules. Indicating with tf^ the hopping am- 
plitude between sites i and j and bands a and/3 and with 
ef the local energy for an electron in the band a on site 
i, our tight-binding Hamiltonian reads 

h = J2 e ? c L w + E ^(4.^ + M+E - 

iacr ijaft i 

(i) 

where Cj a(T and q are creation and annihilation op- 
erators for electrons of spin a on site i and orbital a, 
and rii — J2aa c \aa c ia<r ls the total density on site i. In 
this preliminary study we limit ourselves to onsite (in- 
tramolecular) interactions and we also do not consider 
exchange interactions between LUMO and LUMO+1. U 
controls both interband and intraband interactions. In 
order to implement our mean-field strategy we decouple 
the interaction term in the particle-hole channel, con- 
sidering all possible instabilities (charge and magnetic), 
finding that an antiferromagnetic (AFM) solution with 
order parameter m — Y^ii^iA ~ ^fe) being the z 
component of the spin on the two sublattices) is clearly 
favored for our estimated values of U. Even for the small- 
est estimate given above (U=1.2 eV), the system becomes 
AFM (the A molecules having opposite magnetization 
with respect to B molecules) with a large magnetic mo- 
ment of 0.96 \xb corresponding to (almost) a single spin 
per picene molecule. We notice in passing that a ten- 
dency towards a charge-ordered solution has been ob- 
served within this mean-field approach, even if this solu- 
tion is never stable. Such a tendency is intriguing because 
charge-ordering may be expected in 3/4- filled systems in 
the presence of nearest-neighbor interactions, which are 
not included here. Non-local interactions may stabilize 
this phase, reviving the scenario in which the four bands 
are democratically occupied by the six electrons leading 



to 3/4-filling. 



C. Antiferromagnetic Solution with hybrid 
functionals 

At a mean-field level, we stabilized an AFM state sim- 
ilar to the parent compound of copper oxide supercon- 
ductors or to the ground state of CsaCgo- One can sur- 
mise that also in if-picene superconductivity may appear 
close or compete with an AFM state. It is however well 
known that the mean-field approach emphasizes the ten- 
dency towards antiferromagnetism. Therefore we tested 
this result using HSE and B3LYP hybrid functionals, 
also in light of their success in determining the gap of 
the parent compound that we discussed above. Further- 
more, these approach does not require an estimate of U , 
and represent an independent test of the importance of 
electron-electron correlations. In Fig. [3] we show results 
for HSE, which is believed to perform better in metal- 
lic systems because the long-range tail of the Coulomb 
kernel is screened 2 ^—, even if the results obtained with 
B3LYP are fully compatible. Using both approaches an 
AFM solution with a reduced magnetic moment of 0.4 
Hb at each molecular site is indeed stabilized. Within 
HSE the AFM state has a finite gap of 0.4 eV (see Fig. 
|3p), and the gain in energy respect to the non magnetic 
state is 0.16 eV per unit cell which is at the limit of the 
accuracy of our calculation. B3LYP produces an AFM 
state with the same magnetic moment and slightly larger 
energy gain (0.22 eV) and gap (0.58 eV). 

Therefore, using the same approaches that correctly 
reproduce the gap of the parent undoped compound, the 
trivalent doped system appears to be a low-spin AFM, 
in which less than one electron per picene molecule con- 
tributes to magnetism. This result seems to be close to 
what we obtained by means of static mean-field, with a 
half-filled conduction band, and the smaller value of the 
magnetization can be ascribed either to an overestimated 
U in our mean-field calculation, or to the effect of quan- 
tum fluctuations. Ab-initio estimates of U (and possibly 
of nearest neighbor repulsion and interorbital exchange) 
will be very helpful to understand this discrepancy, as 
well as more accurate treatments of short-ranged corre- 
lations, such as Dynamical Mean-Field Theory. 



IV. CONCLUSIONS 

In this manuscript, motivated by experiments showing 
that for 2.7 < x < 3.3 if^-picene displays superconduc- 
tivity with T c = 18K, we have investigated the band- 
structure of this compound with a special focus on the ef- 
fects of electron-electron correlations. In agreement with 
previous DFT calculations, we find that the potassium s- 
electrons are donated to four bands arising from LUMO 
and LUMO+1 of two inequivalent picene molecules. Six 
electrons are donated in the case of i^-picene. This gives 
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rise, within our DFT methods, to four completely filled 
bands arising from the LUMOs while the two bands aris- 
ing from the LUMO+ls are overall half-filled. Estimates 
of the correlation strength pose these materials on the 
correlated side, suggesting that magnetic ordering can 
be expected. While the tendency toward magnetism is 
very weak in PBE, the inclusion of correlations at the 
mean-field level gives a strong tendency towards antifer- 
romagnetism, which is confirmed by the accurate hybrid 
functional HSE (and also by B3LYP), even if the mag- 
netic moment is substantially reduced. The antifcrro- 
magnetic ordering seems quite robust, and it should not 
depend on the details of the crystal structure, which is 
not known experimentally for the doped system. 

This finding would put -fCr-picene in a wide class of 
superconductors in which AFM and superconductivity 
coexist and/or compete in the same phase diagram. On 
the other hand, a very recent photoemission study2£ sug- 
gests that a sizable intramolecular electron-phonon cou- 
pling is present in i^-picene with x ~ 1, and a role of the 
same kind of phonons has been proposed in a recent DFT 
analysis^. This is particularly interesting in light of the 
analogy with fullerides. In these compounds an intra- 
molecular electron-phonon coupling is indeed believed to 
provide pairing^, and it has been shown that a simi- 
lar coupling can indeed be favored by strong Coulomb 
interaction 4 - because it does not involve charge fluctua- 
tions on a molecule, but rather it couples to the internal 



degrees of freedom, which are not frozen by Coulomb 
interactions. It appears therefore that this new molecu- 
lar superconductors shares many similarities with alkali- 
metal doped fullerides. In this light the AFM state that 
we find in the present study would be competing with 
a phonon-driven (and correlation assisted) groundstate, 
leading to a first-order transition as a function of the cell 
volume, in close analogy to what observed in CsaCecn. 

After completion of the manuscript, we became aware 
of a related work by Kim et at, arXiv:1011.2712 ( in which 
an AFM state is found also using PBE if the volume per 
molecule is increased by 5%. This result is compatible 
with ours, and confirms the important role of electron- 
electron correlations in i^-picene. 
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